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Abstract 

In the kinematical region where the center of mass energy is much larger than ah other 
scales, the Drell-Yan process can be formulated in the target rest frame in terms of 
the same color dipole cross section as low Bjorken-x deep inelastic scattering. Since 
the mechanisms for heavy dilepton production appear very different in the dipole ap- 
proach and in the conventional parton model, one may wonder whether these two 
formulations really represent the same physics. We perform a comparison of numerical 
calculations in the color dipole approach with calculations in the next-to-leading order 
parton model. For proton-proton scattering, the results are very similar at low X2 
from fixed target to RHIC energies, confirming the close connection between these two 
very different approaches. We also compare the transverse momentum distributions 
of Drell-Yan dileptons predicted in both formulations. The range of applicability of 
the dipole formulation and the impact of future Drell-Yan data from RHIC for deter- 
mining the color dipole cross section are discussed. A detailed derivation of the dipole 
formulation of the Drell-Yan process is also included. 
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1 Introduction 



With the advent of RHIC, the Drell-Yan(DY) process []T| can be studied in a new kinematical 
regime, the so-called Regge regime, where the dilepton mass M is small compared to the cm 
energy -y/i, but still much larger than Aqcd- The DY process at RHIC and LHC energies is 
therefore of similar interest as DIS at HERA, where one can study 7*-proton scattering in the 
Regge regime. The new experimental possibilities motivate further theoretical investigations 
of the DY process. First of all, one needs a framework to calculate nuclear shadowing in the 
DY process, the onset of which can already be observed at fixed-target energy P|. This is 
especially important in view of the RHIC heavy ion program. The color dipole formulation 
of the DY process introduced in j^, Q is suitable to address this issue and a considerable 
amount of work in this direction has been performed 0-0. In addition, the I0W-X2 DY cross 
sections are sensitive to integrals over the color dipole cross section which are not accessible 
in DIS [Q. Therefore, future DY data can be used to further constrain this quantity. Most 
importantly, while the dipole approach has already been used in sophisticated analyses like 
the extraction of energy loss from E866/NuSea Drell-Yan data f^, its validity for proton- 
proton interactions has never been established. This will be remedied with the present work. 

In this paper we compare next-to-leading order (NLO) parton model calculations for DY 
dilepton production in proton-proton (pp) and proton-deuteron (pd) collisions with calcu- 
lations in the dipole approach over a wide energy range. Although the two approaches are 
believed to be equivalent in a certain kinematical range, the underlying mechanisms appear 
to be quite different, and there is no known way to prove this equivalence analytically. How- 
ever, both approaches are supposed to describe the same process, so they should yield similar 
numerical results^. 

Before we compare the results of the numerical calculations, the key features of the two 
approaches are briefly summarized. The well known mechanism for continuum dilepton 
production, which was first found more than thirty years ago by Drell and Yan was 
formulated in a frame where both colliding hadrons are fast moving (infinite momentum 
frame). According to Feynman's picture of high energy collisions, the colliding objects can 
be viewed as collections of noninteracting partons with negligibly small transverse momenta. 
To lowest order, DY dileptons are produced by quark-antiquark annihilation, and the cross 
section reads, 

Nf 

The distribution function of a quark (antiquark) of flavor / of the target or the projectile 
is denoted by qf {qf), Nf is the number of active flavors and Zf is the quark charge. The 

^Differences between parton model and dipole approach are seen in the angular distribution of DY pairs, 
if saturation effects are included in the dipole cross section ||l0| . When we talk about equivalence of parton 
model and dipole approach, we mean equivalence up to higher twist effects, which are neglected in the parton 
model, but (at least partially) included in the phenomenological parameterization of the dipole cross section 
we employ. 
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longitudinal momentum fractions of the projectile (target) parton, Xi (X2), can be expressed 
in terms of Lorentz invariant scalar products as 



2P2-g 2P,-q 

xi = ; X2 = , (2) 

s s 

where {P2) is the projectile (target) four momentum, is the four momentum of the 
dilepton, = > 0, and xp is the Feynman-x, xp = xi — X2- 

For most qualitative descriptions, it is sufficient to consider the DY process in terms of 
the lowest order annihilation process, Eq. (|I]). Calculations with Eq. (|l|), however, under- 
estimate measured DY cross sections by an overall factor. It is necessary to employ the 
NLO framework for the DY process, in order to make quantitative predictions, see for a 
review. In addition, the DY cross section differential in the dileptons transverse momentum 
receives huge corrections from higher order processes. Indeed, to lowest order, one would 
not expect dileptons with large transverse momentum q±, in contrast to what is observed 
in experiment. Even though the occurrence of perturbatively large transverse momenta can 
be explained in NLO, it is not straightforward to calculate the shape of the g^-distribution 
in the parton model. A resummation of large logarithms in q±/M |jl2| or alternatively the 



introduction of an intrinsic transverse momentum |T3| is necessary to avoid the divergence 
of the differential cross section at g_L = 0. 

In the parton model, all nonperturbative effects are parameterized in the parton distri- 
bution functions qj, qj, which evolve according to the DGLAP evolution equations. For 
DY in nuclear collisions, the parton distribution functions of the proton are simply replaced 



by empirical nuclear parton distribution functions ||14|. This approach does not explain the 
dynamical origin of the nuclear effects 




Figure 1: In the target rest frame, DY dilepton production looks like 
hremsstrahlung . A quark or an antiquark from the projectile hadron scat- 
ters off the target color field ( denoted by the shaded circles ) and radiates 
a massive photon, which subsequently decays into the lepton pair. The 
photon decay is not shown. The photon can be radiated before or after 
the quark (antiquark) scatters. 

Nuclear effects, effects from higher orders in perturbation theory, as well as other possible 
nonperturbative effects, are more readily treated when the Drell-Yan process is viewed in 
the target rest frame. Note that although cross sections are Lorentz invariant, the partonic 
interpretation of high energy scattering processes does depend on the reference frame. In the 
rest frame of the target, the production mechanism for high mass continuum dileptons looks 
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like bremsstrahlung see Fig. |l|. In the high energy hmit, when one can neglect terms 

that are suppressed by a factor 1 / energy, each of the two graphs factorizes into a production 
vertex for the virtual photon times an amplitude for scattering a quark off the target. These 
scattering amplitudes combine in the squared matrix element in exactly the same way as 
in DIS, which makes it possible to express the DY cross section in terms of the same cross 
section cr^ for scattering a gg-dipole off a nucleon (A^) as in low-x^j DIS, 



din a 



d p \^^,q{a,p)\ agg{ap,x) . (3) 



Here, a is the light-cone momentum fraction the virtual photon takes away from its parent 
quark, and p is the transverse separation between 7* and final quark. The electromagnetic 
radiation, q -h> 7*g, is described by the light-cone wavefunction \l/^*g(a,p), see Eqs. ( |A.18| ) 



- ( |A.20| ), which can be calculated perturbatively. Summation over photon polarizations is 
understood in Eq. (|^). The dipole cross section cx^ is of nonperturbative origin and has to be 
taken from phenomenology. The energy scale x of the dipole cross section will be discussed 
in the next section. A detailed derivation of Eq. @ is given in the appendix. 

Using a phenomenological parameterization for the dipole cross section in Eq. (|^) is a 
very economical way to account for higher order and nonperturbative effects. The dipole 
approach can even be applied at low values of M where perturbative QCD is not valid [|I5|. It 
was found in a recent analysis [|10| that most of E772 DY data (except some points at low M) 



are reasonably well described in the dipole approach without introducing an arbitrary overall 
normalization factor. In addition it was found that the transverse momentum distribution 
does not diverge at q± = 0, even without intrinsic transverse momentum. 

We emphasize that the dipole approach does not describe an additional production mech- 
anism for heavy dileptons. Rather, the two approaches are believed to describe the same 
physics in different reference frames. Therefore, calculations in the NLO parton model and 
in the dipole approach should give similar results for the DY cross section. This is what we 
numerically check in this paper. In the following section, we compare numerical calculations 
of the DY cross section (integrated over the transverse momentum of the dilepton) in both 
approaches. In section 3, we also compare the predictions of dipole approach and parton 
model for the DY transverse momentum distribution at RHIC. 



2 Numerical comparison of the two approaches 

In order to perform calculations that can be compared with experimental data, one has 
to embed the partonic cross section, Eq. (^), into the hadronic environment. In the infinite 
momentum frame, the momentum fraction of the projectile quark is Xi, see Eq. (|]). However, 
when the scalar product defining xi is evaluated in the target rest frame, one finds Xi = az, 
where z = Xi/a is the momentum fraction of the incoming proton carried by the projectile 
quark. The different meanings of xi in the target rest frame and in the infinite momentum 
frame is a manifestation of the frame dependence of partonic mechanisms. In the target rest 
frame, xi is the momentum fraction that the lepton pair takes from the projectile proton. 
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Thus, one obtains for the proton-nucleon DY cross section 



(Pa{pN I'^l X) OL(,m Xi da 



1 f=l 



dW^dxF 37rM2 xi + X2 

X j d^p\m^,^{a,p)f a^q{ap,x) . (4) 

We still need to know the scale Q at which the projectile parton distributions are probed 
and the energy x at which the dipole cross section enters. These quantities are not known 
exactly, instead we have to rely on plausible arguments to determine their values. In order 
to find Q, note that the transverse distances p that contribute to the DY cross section are 
controlled by the extension parameter 

rf = {l-a)M'^ + a^m). (5) 

The numerically dominant term in the LC wavefunctions, Eqs. ( |A.18| , |A.19| ), is the one 
that contains the Bessel function Ki(r7p). Since this function decays exponentially at large 
arguments, the largest distances that can numerically contribute are of order ~ I/?]. For 
fiuctuations with a — > 1, these distances can become of the order of a typical hadronic radius, 
in analogy to the aligned jet configurations in DIS On the other hand, the minimal value 



of a is Xi, so that the largest virtuality entering the calculation is = ry^ax = (1 ~ Xi)M^. 
We choose this quantity to be the hard scale at which the projectile parton distribution 
is probed. The parton distribution functions (PDFs) are taken from CERNLIB |]T^. The 
quark mass is set to = in all our calculations, see for its numerical infiuence. 

For the quark density of the projectile, we employ the leading order parameterization that 
corresponds to the NLO parameterization used in the parton model calculation. This means 
e.g. we use CTEQ5L in the dipole approach when comparing it with a NLO parton model 
calculation using CTEQ5M. One should use leading order PDFs in the dipole approach, 
because they are scheme independent and have a probabilistic interpretation. 

The energy scale x of the dipole cross section in Eq. (^) is determined from the analogy 
to DIS. In DIS, the argument of the dipole cross section is xbj = Q'^/W'^, where Q is 
the virtuality of the photon and W is the 7*-proton cm energy. Therefore, we choose 
X = M'^/s = ax2, where s = sxi/a is the quark-proton cm energy squared. 

Note that in the previous analysis [jlO|, and X2, instead of and x, were used. The 
different choice of scales in this paper has the effect of increasing the cross section by a factor 
of up to 2 for dilepton mass M ~ 4 GeV. This is mostly due to the different choice of Q^. 
Using ax2 instead of X2 is only a ~ 10% effect at X2 < 0.1. These uncertainties vanish at 
larger masses, M ~ 8 GeV. 
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Figure 2: Calculations in the dipole approach to DY and its modifica- 
tion Eq. (|^ compared to NLO parton model results at fixed target energy 
(y/s = 38.8 GeV). The CTEQ5M parameterization J^/ is used for the 
parton model calculation. The data are from E772 The curves and 
data for the different mass bins have been rescaled for better visibility. In 
all calculations, none of the free parameters of the dipole approach were 
adjusted to DY data. Only DIS data have been used. 
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For the calculations in the dipole approach shown in Figs. - we employ the parame- 
terization of the color dipole cross section by Golec-Biernat and Wiisthoff p 



1 — exp 



(6) 



with (To = 23.03 mb. This parameterization rises quadratically at small separations p as 



demanded by color transparency |T9[ and saturates at large separations. The saturation 



scale that controls the flattening of the dipole cross section is given by 

We point out that no DY data have been used to determine the parameters in Eq. (P). Only 
DIS data from HERA were fltted to extract the dipole cross section. 

Note that the parameterization Eq. (|^) contains only the Pomeron part of the dipole 
cross section. As a consequence, the dipole approach predicts identical cross sections for 
particle and antiparticle induced DY. This is, of course, reproduced in the parton model 
at low X2- In addition to this Pomeron part, there is also a contribution from Reggeon 
exchange, corresponding to valence quarks in the target |^. In principle, one could flt 
a parameterization of this Reggeon part to the data, but the predictive power of such a 
procedure would be very low, since the Reggeon part, in contrast to the Pomeron part, 
depends on the colliding hadrons. We therefore do not attempt to determine the Reggeon 
contribution to the dipole cross section in this paper. 

The NLO parton model calculation was performed with the CTEQ code pO| provided 
earlier to the E866/Nusea collaboration by W.K. Tung [^. We use the dilepton mass M 



as hard scale in the strong couphng constant and in the parton distributions. This is the 
standard choice. The code assumes MS renormalization scheme. 

Fig. shows calculations at flxed target energy in the dipole approach and in the NLO 
parton model, as well as experimental DY data from E772 [^. The lowest value of X2 
that can be reached is about 0.02. Dipole and parton model calculations are very similar at 
low X2, where the dipole approach is supposed to be valid, but strongly deviate as X2 — 1. 
Clearly, we cannot apply the dipole approach at large X2, since it does not include valence 
quark contributions to the DY cross section (see |jlO[ for a more detailed discussion of effects 
that are taken into account by the dipole approach). Note also that the parameterization 
Eq. was fltted only to DIS data below xbj < 0.01, hence the curves shown in Fig. ^ are 
already an extrapolation of the flt. 

In Fig. there are several E772 data points at low X2 that exceed the calculations in 
both the parton model and dipole approaches. This excess was also visible for some points 
in Indeed, experiment E772 was not designed to measure the absolute normalization of 
the DY cross section, but to measure nuclear effects. New data from E866/NuSea will not 
suffer from such problems. A comparison between preliminary E866/NuSea data and the 



NLO parton model in shows good agreement 



The disagreement between the two approaches sets in around X2 ~ 0.1 in Fig. ^. Some- 
what surprisingly, the dipole approach yields values that exceed the prediction from the 
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parton model (at large X2), even though the dipole approach does not include several con- 
tributions to the DY cross section. This can be understood, if one observes that the dipole 
cross section is related to the target gluon density G. At small separations p one has 



<(p,x) = y«.(A^AG(x,A). (8) 



Eqs. (P) and (^) contain only the part of the gluon density that increases at low x. At large 



X, however, the gluon density decreases like G{x) oc (1 — x)^ At high virtuality, the 

power is even higher than 5, say ~ 8. In order to estimate the uncertainty originating from 
this effect, we follow and replace the saturation scale Eq. (0) by 



Qiix)-^Qi{x){l-xr. (9) 

We refer to this recipe as modified dipole approach and the numerical results are shown by 
the short dashed curves in Figs. ^ The replacement Eq. (^) has virtually no influence 
at X < 0.01, where the parameters in Eq. (^) have been fitted in [|18|, but at x = 0.1, it 



reduces the cross section by almost a factor of 2. In fact, in [^] calculations in the dipole 
approach were performed employing a dipole cross section calculated from a gluon density 
that contains saturation effects. For y/s = 38.8 GeV, the results of ||2^ lie below the ones 



obtained in |T0[ with the saturation model of Golec-Biernat and Wiisthoff [|T^ by a factor of 
~ 2. We believe this disagreement is largely due to the decreasing gluon density at X2 > 0.01, 
i.e. due to uncertainties in the dipole cross section. 

Apparently, the I0W-X2 range accessible at fixed target energies is very limited, and there- 
fore the dipole formulation is bestowed with several uncertainties. Unfortunately, there are 



no DY data at really low X2- Even CDF data [|30| are at about the same values of X2 as 
the E772 data, because they are mostly in the Z-boson mass range. At RHIC, however, 
much lower values of X2 will be reached, where the dipole approach can readily be applied. 
In Fig. ^ we show predictions from the dipole approach and the NLO parton model (using 
CTEQ5M PDFs) for RHIC energy. Calculations with the modified dipole cross section, 
Eq. (P), are also shown. The disagreement between dipole approach and parton model sets 
in around X2 ~ 0.01, i.e. around xp ~ 0. At low X2 all three curves are in good agreement, 
confirming the similarity between the two approaches. 

In order to show the uncertainty arising from the choice of PDFs, we did the same 



calculation, employing GRV98HO [|2J], see Fig. ^. In this case, the differences between the 
two approaches are slightly larger than for CTEQ5M. With GRV98HO, the NLO parton 
model calculation exceeds the dipole approach by at most a factor of 1.5 for X2 < 0.01; for 
CTEQ5M, the corresponding value is only about 1.2. The third standard PDF we employed 
was MRST(c-g) ||31]|. The results are comparable to GRV98HO and are not shown here, 
because they do not provide any further insight. For most values of X2 and M, which do 
not lie at the very edge of the phase space, the NLO calculations are only ~ 10% - 20% 
higher than the dipole approach. In summary, we find that the uncertainty from the choice 
of PDFs is ~ 25% at the low values of X2 shown in Figs. ^ and |^. 
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Figure 3: The same as Fig. but now at RHIC energy (^/s 
GeV), where much lower values of X2 can be reached. 
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Figure 4: The same as Fig. ^ but with GRV PDFs j^/ . For the dipole 
approach, we employed the LO PDFs GRV98LO for the projectile parton 
distribution. 
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We also point out that for these low values of X2, the DY i^-factor is typically 1.4 to 1.5 
and can become even smaller than unity for very low X2- In the parton model, the i^'-factor 
is defined as the DY cross section calculated in NLO divided by the same quantity calculated 
in LO. The /^-factor originates (mostly) from the analytic continuation from the spacelike 
(DIS) to the timelike (DY) region of in the annihilation and virtual corrections at NLO. 
The QCD Compton process, which dominates at low X2 because of the large target gluon 
density, does (almost) not contribute to the J'T-factor. Therefore, in the kinematical region 
of interest, the i^-factor is considerably below its usually assumed value of 2 to 3. 

At this point, we stress that the term "i^-factor" is defined only in the parton model. 
In the dipole approach, there is no freedom to adjust the overall normalization of the DY 
cross section as proposed in [^, because higher order corrections are contained in the 
parameterization of the dipole cross section. 

Most uncertainties in the dipole approach arise from uncertainties in this parameteriza- 
tion. It remains to be seen whether future low X2 DY data from RHIC can be accurately 
described with a (probably improved) phenomenological parameterization of the dipole cross 
section. We believe that future DY data can serve as an important source of knowledge about 
the color dipole cross section. 



3 The DY transverse momentum distribution 

The DY cross section differential in the dilepton transverse momentum q± can provide even 
more detailed information about the shape of the dipole cross section than the g_L-integrated 
cross section. In the dipole approach, the g^-differential cross section reads P], 



d?a{qp 7*X) 



d\n.ad?q± (27r 



^— j dViC?V2e^^^-('^^-^^^)^;.,(a,Pi)vl/,., 



X ]^{aqq{api,x) + aqq{ap2,x) - aqq{a{pi- p2),x)] , (10) 

where q±_ is the 7* transverse momentum wrt the direction of the projectile quark. This 
formula was first published in 0. We give a detailed derivation in the appendix. The 
partonic cross section, Eq. (|l^), has to be embedded into the hadronic environment in the 
same way as in Eq. (^). Again, the sum over photon polarizations is understood to be 
contained in the light-cone wavef unctions, see Eq. ( |A.20| ). 

We would like to compare the transverse momentum distribution calculated in the dipole 
approach with the same quantity computed in the parton model. Since the CTEQ code 
used in the previous section is not capable of calculating the DY transverse momentum 



distribution, we developed our own code based on the formulas given in [0. Again, we use 
the dilepton mass M as hard scale. It was already mentioned in the introduction, that there 
are no large transverse momenta in the leading order parton model. The order as correction 
yields a transverse momentum distribution that diverges at g_L ^ and has the wrong shape, 
when compared to data |TB|. This can be remedied by resumming contributions from soft 
gluon radiation |12| , which yields a good description of the data. 
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In this paper, however, we follow |13| and apply a more phenomenological recipe. We in- 
troduce a soft, nonperturbative, primordial transverse momentum distribution of the partons 
in the colliding protons, which we parameterize by a Gaussian, 



The mean transverse primordial momentum squared of a single constituent in the proton is 
then 

primordial = 2(7^ . (12) 

The perturbatively calculated g^-distribution is then smeared out by this primordial mo- 
mentum and the regularized dilepton cross section is given by 

/ f]2 DY \ 
\dM'dXF J ^^O+LO 

The subscript A + C refers to "annihilation-|-Compton" correction. The virtual corrections 
do not contribute to the first term in Eq. (p!3[). With this recipe, the perturbative result is 
reproduced for large transverse momenta while the divergence at g_L = is removed. The first 
term on the rhs of Eq. (|T3|) vanished after integration over q± and J d?q±f{q±) is normalized 
to unity, so that the total DY cross section is correctly reproduced. 

The results of our calculations in the dipole approach and in the parton model are shown 
in Fig. ^. We performed the parton model calculation for two different values of the smearing 
parameter cXg. The lower value, Ug = 0.48 GeV, was used in |^ to describe ISR data at 
^/s = 44 GeV and 62 GeV. The parameter aq is presumably independent of energy s. 
Nevertheless, we also calculated for twice the ISR value of aq, because at the low values of x 
reached at RHIC, quarks may have a larger intrinsic transverse momentum. Of course, the 
choice of cxg has the largest influence at small q±, while at very large transverse momentum, 
the results become independent of this parameter. Note however, that aq still has a rather 
large impact on the numerical results at intermediate ~ 2 - 3 GeV. 

No intrinsic transverse momentum is included in the dipole approach, even though one 
could introduce one. As already pointed out in 0], the cross section calculated in the dipole 
approach does not diverge at — > because of the flattening of the dipole cross section, see 
Eq. d^). At q± = 10 GeV, dipole approach and parton model differ by about a factor of 5. 
Employing the modified dipole approach, Eq. (§), brings the curve down by a factor of 0.9 
at q± = 10 GeV and has no effect at small q±. Therefore, we do not show the corresponding 
curves. Note that at high transverse momentum, the numerical value of the DY cross section 
is completely insensitive to the large p behavior of the dipole cross section. Indeed, it was 
found in ||10[ that an expansion of Eq. to order virtually yields the same numerical 
results at large q± when employed in Eq. (plQl). On the other hand, the behavior at 
is largely determined by the large p behavior of the dipole cross section. Therefore, future 
data on the DY g^-distribution can be important for further constraining the dipole cross 
section, provided statistics will be high enough. 
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Figure 5: Comparison of the DY transverse momentum distribution 
calculated in the dipole approach and in the NLO parton model. The 
parton model calculations have been performed for two different values of 
the quark's intrinsic momentum. The values of the kinematical variables 
for all curves in this figure are ^/s = 500 GeV, xp = and M = 4.25 
GeV. 

4 Summary 

In this paper, we presented a detailed comparison between calculations of the Drell-Yan cross 
section in the color dipole approach with calculations in the conventional NLO parton model. 
We also compared both approaches with experimental data when available. The dipole 
approach and the parton model are believed to represent the same physics, but viewed from 
different reference frames. The purpose of this paper was to demonstrate the close connection 
between the two approaches numerically. 

First, we compared calculations for the DY cross section integrated over the dilepton's 
transverse momentum q± and found that at low X2, where the dipole approach is supposed 
to be valid, calculations agree well. Since both approaches rely on different nonperturbative 
inputs, numerical results will not agree exactly. There are several uncertainties from free 
parameters that can affect the numerical value of the Drell-Yan cross section. In the dipole 
approach, these uncertainties are: 

• Different possible choices of the hard scale Q, at which the parton distribution of the 
projecetile is probed, can lead to an uncertainty as large as a factor of 2 for M ~ 4 GeV, 
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even though the dependence on Q is only logarithmic. This uncertainty vanishes with 
increasing dilepton mass M . For M = 8.75 GeV it is only ~ 10%. Note that a larger 
value of Q leads to a smaller value of the DY cross section, because the projectile PDF 
is probed at large momentum fraction. 

For the low values of X2 we are interested in, different parameterizations of PDFs lead 
to variations of the numerical results of up to ~ 25%. This uncertainty is also present 
in the parton model. 

The uncertainty arising from possible choices of the energy scale x is ~ 10%. Since the 
dipole cross section increases with decreasing x, a lower value of x pulls the DY cross 
section up. 

The dipole approach allows one to introduce a quark mass ruf. In this work we set 
rriq = 0. A quark mass = 200 MeV would reduce the numerical value of the DY 
cross section by less than 10% ITU . 



The main uncertainty, which is difficult to quantify, arises from the parameterization of the 
dipole cross section. Future DY data from RHIC will certainly help to better constrain this 
quantity. 

We found that the disagreement between the two approaches is typically of order 20% 
for X2 < 0.01. We believe that the parameterization of the dipole cross section can be 
adjusted so that future RHIC data can be described without introducing an arbitrary overall 
normalization factor for which we cannot find a theoretical justification. For larger values of 
X2, the applicability of the dipole approach is questionable, because this approach neglects 
several contribution to the DY cross section which might already be important for 0.01 < 

X2 < 0.1. 

Note that several data points from E772 at low X2 are not well described by neither 
of the approaches. Preliminary E866/NuSea data p5[, however, which were measured in 



the same kinematical region, are well described by the NLO parton model. Therefore, the 
agreement between dipole approach and NLO parton model at low X2 gives confidence into 
the applicability of the dipole approach. However, we stress that in order to strictly prove 
the validity of the dipole approach, one has to reproduce experimental data and not the 
parton model calculation. 

Naturally, the differences between the two approaches are larger when one considers 
the DY transverse momentum distribution. At low q±, the numerical result in the parton 
model depends quite strongly on the amount of primordial transverse momentum which is 
included. This uncertainty is however no longer present at very large transverse momenta, 
gj_ ~ 10 GeV, where the two calculations differ by a factor of 5 (at y/s = 500 GeV) . It 
will be seen, if RHIC data can distinguish between the two approaches or if the dipole cross 
section and/or the parton distribution function can both be adjusted to fit the data. We 
point out that both approaches predict the same asymptotic behavior for the partonic DY 
cross section, namely da/d'^q± oc 1/g^ for large transverse momentum 
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Appendix A Derivation of the dipole formula for DY 

In this appendix, we give a detailed derivation of the formula that expresses the DY cross 
section in terms of the color dipole cross section, Eq. (|]). This formula was first discovered 
in 0]. Some details of the derivation in the semiclassical approach of |Q are given in the 
appendix of 0]. In this paper, we use the framework of covariant Feynman perturbation 
theory. The calculation can be carried out in a very similar way for DIS. 

The differential cross section for dilepton production in the target rest frame, see Fig. ^ 
is given by 

xxx^ (27r)58(p^)2a(l - a) 

X „.„..(A).;(A')L-i^, (A.1) 

where pi is the momentum of the projectile quark, p/ is the momentum of the final quark 
and q is the momentum of the massive photon, q^ = > 0. The first line describes the 
bremsstrahlung of a massive photon from a quark (as depicted in Fig. [ip, and the second 
line the decay into the lepton pair into solid angle dQ (not shown in Fig. 0). The 
electromagnetic coupling constant is denoted by Oem = e^/ (47r) = 1/137. Furthermore, e^(A) 
is the polarization vector of the photon, A G {±1,0}. We sum over all final states X and 
the photon polarizations A, A'. The 5-function for energy conservation is already integrated 
over. Furthermore, a is the energy fraction the photon takes from its parent quark. It is 
assumed that only transverse momentum but no energy is exchanged with the target. Then, 
Ai^'^ is the absolute square of the bremsstrahlungs-amplitude, before contraction with the 
polarization vector of the photon. There are two contributions to this amplitude, 

^'^ = IJ2w^ + -^umT + Min. (a.2) 

afai CfCi 

The s-channel amplitude describes the process where the photon is radiated after the 
quark has scattered off the target. Fig. |I| (left), and the w-channel amplitude describes 
the process where the photon is radiated before the quark scatters off the proton. Fig. |l] 
(right). The bar indicates that the matrix element is summed over helicity aj and color Cf 
of the final quark and averaged over helicity cTj and color Cj of the initial quark. (A*",, = 3 is 
the number of colors.) Finally, pi± are the dilepton momenta and 



' 111/ 



^{p^^Pi- + Pi^P^- - g'^Pim- ) ( A.3) 
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is the leptonic tensor, summed over helicities. 

Different terms in the sum over photon hehcities lead to different angular distributions of 
the lepton pair. In this paper, we integrate over the solid angle dfl = d(pd{cos 9) , and obtain 
for the differential DY cross section 



d\iiadM'^d'^qj_ SvrM^ |^ d\iiiad'^q± d\nad'^q± 

where the indices T and L stand for transverse and longitudinal photons, respectively. We 
also integrated over the phase space of the final quark. The factor aem/(3vrM^) originates 
from the second line of Eq. ( [A.l| ) and describes the decay of the 7* into the lepton pair. The 



interesting physics resides in the quantities d^aT^i/dlna d'^q±, which are given by 
d^ariqN^rX) _ f ^ ^ 6* (A)e,(A)7W' 



X AG{±1} 



dlnad-'q^ J ^ {27if8{p°,y{l-a 



Terms with A 7^ A' in Eq. ( [A. 1|) vanish after integration over the azimuthal angle 



We will now express the d^aT.b/ d\\iad'^q± in terms of the q —>■ '~f*q light-cone (LC) wave- 
functions and the dipole cross section. Consider the s-channel graph in Fig. |l] (left) where 
the photon is radiated after the projectile quark has scattered off the target. The propagator 
of the quark in the intermediate state can be written as 

^f+i + mj _ ^ Ua{pf + q)ua{pf + q) 7 



{pf + qy — m'j ^ [pf + qy — m'j 2(pjr + g+)' 

The u^{p) are Dirac-spinors, {^—mf)u„{p) = 0, for on-shell momentum p and helicity a. The 
7"^-term arises, because the intermediate quark is off-shell. In the high energy approximation 
which we employ, this term is dropped^. This is the crucial step that allows one to write the 
amplitude as a product of an LC wavefunction and a quark scattering amplitude. 
The s-channel amplitude then turns out to be 

For simplicity, we have set the flavor charge to unity. Here, tg,o-o-i((P/ + q^),k±) is the 
amplitude for scattering a quark off a nucleon in the rest frame of the nucleon, 

t,,..Aip'f + q'), k±) = u,{pf + q) 7° V,{k^) u^M - 2pX,^y,{k^). (A.8) 



■^As already pointed out in [Q, this approximation is equivalent to neglecting the instantaneous vertices 
in a Hamiltonian light-cone approach to QCD | |3^ , and violates gauge invariance. The color dipole approach 
is gauge invariant only in leading log(a:;) approximation. 
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Note that t and V are matrices in color space, we suppressed the color indices. Our spinors 
are normalized to u^^{p)u^i{p) = 2p°(5o-,o-'- The function Vg{k±) will be completely absorbed 
into the dipole cross section, so that we will never need to specify it. Therefore, the dipole 
formulation also accounts for some higher order and nonperturbative effects, which are dif- 
ficult or even impossible to account for in the standard parton model. We write Vq only as 
function of the exchanged transverse momentum k± = Pf± + q± — Pi±, because the longi- 
tudinal momentum of the projectile quark cannot be changed significantly at high energies. 
Note that Vq also depends on energy, even though we do not write this dependence explicitly. 
The Kronecker-(5o-,(7i means that the helicity of the quark is not changed by scattering off the 
proton. 

In impact parameter space, the amplitude reads 

M^sib^p) = j '^^^^^-''^■'-'^''^^^^ (A.9) 

= -iv/4;^^^lH^v[;/; («,p)2pOr,(6), (A.IO) 

where l± = Pf± — {l — (y)q±/<y is the transverse momentum of the final quark in a frame where 
the 2;-axis is parallel to the photon momentum. The conjugate variable to l±, namely ap, 
is then the transverse distance between initial and final quark in that frame. This quantity 
will become the distance variable the dipole cross section. Furthermore, 

VM = J (A.ll) 

The LC wavefunctions in impact parameter space are related to the quark-photon vertex 
and the propagator by Fourier transformation, 

,TrM f ^ f d^^^ ~il^-ap , U .{pf)-f''u^XPf + Q) . . ^ 

%,q{a,p) = aWl-a J ^ aHl + r]'' ' 

where 7f = {1 — a)M'^ + a'^m'j. In our numerical calculations in this paper, we set the quark 
mass ruf = 0. See [0] for its influence. 

The M-channel graph can be written in the same way as ( |A.9| ), but with a quark scattering 
amplitude at a shifted impact parameter. This can be seen from the propagators. While the 
propagator for the s-channel graph yields 

^ "^^-"^ (A.13) 



one finds for the u-channel graph the combination l± + k±_ instead of 

1 a 

{pi - qf -mj~ a\ti_ + k^Y + r/2 ' 



(A.14) 
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One then obtains for the w-channel amphtude in impact parameter space 



= iv/4^^^^vl>i;.^(a,p)2pOy,(6 + ap-). 
and with Eqs. ( [A .41) and ( |A.5| ) for the photon production cross sections 



(A.15) 
(A.16) 



(fpf± 



(2vr) 



(A.17) 



X CfCi 



where l± = pf± — (1 — a) q±/a and /cj, = pf± + q± — Pi±. In the transverse case (T), the 
sum over A includes both transverse polarizations, A = ±1, while in the longitudinal case 
(L) only one term, A = appears in the sum. 

The LC wavefunctions for radiation of transverse and longitudinal photons are readily 
calculated by contracting Eq. (|A.12| ) with the photon polarization vector |, ^, 

^'^.^(a,Pi)^^.;(a,P2) = ^ i^e;(A)^!;.^(a,pi)e^(A)^!;:^(a,P2) 

A=±l CTfCT; 



Op 



em I 2 4 



27r2 



rrija Kq (?7Pi) Kq (?7P2 



[l + (l-af]ry^^^Ki(r/pi)Ki(W2) 1^, 



(A.18) 



C^em ^^2/1 _ \ 2 



M^(l-a)^Ko (r/Pi)Ko (r/p2) 



(A.19) 



Often, one needs only the sum over all polarization states of the photon. We denote the 
corresponding quantity by 

pi)^;*g(a, P2) = Pi)^%{a, P2) + Pi)^^*q{a, P2). (A.20) 

Integrating ( |A.17| ) over the transverse momentum q±_ of the photon yields 



dlna 



d'p 



(A.21) 
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with the dipole cross section 



<?<? 



V,{h)-V,{h + ap) 



(A.22) 



If one assumes that the quark interacts via one-gluon exchange with the target and expands 
Eq. (|A.22|) to first order in (ap)^, keeping only terms of order a^, then the dipole cross 
section is proportional to the target gluon density 

We can also express the transverse momentum distribution of DY pairs in terms of the 
dipole cross section [^. The differential cross section is given by the Fourier integral 



(27r) 



(A.23) 



To derive this expression, one performs the integration over pf±^ in Eq. (|A.17|) and observes 
that Eq. (|A.17|) has a real value. This allows one to symmetrize the integrand with respect 
to pi and p2- The functions Vq combine then to the dipole cross sections in the second line 
of Eq. (|A.23|) . After integrating Eq. ( |A.23|) over the transverse momentum g_|_ of the photon, 
one obviously recovers Eq. ( |A.21| ). 
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